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Abstract 

The physical origin is investigated of Robin boundary conditions for 
wave functions at an infinite reflecting wall. We consider both Schrodinger 
and phase-space quantum mechanics (a.k.a. deformation quantization), 
for this simple example of a contact interaction. A non-relativistic particle 
moving freely on the half-line is treated as moving on the full line in 
the presence of an infinite potential wall, realized as a limit of a Morse 
potential. We show that the wave functions for the Morse states can 
become those for a free particle on the half-line with Robin boundary 
conditions. However, Dirichlet boundary conditions (standard walls) are 
obtained unless a mass-dependent fine tuning (to a reflection resonance) 
is imposed. This phenomenon was already observed for piece- wise flat 
potentials, so it is not removed by smoothing. We argue that it explains 
why standard quantum walls are standard. Next we consider the Wigner 
functions (the symbols of both diagonal and off-diagonal density operator 
elements) of phase-space quantum mechanics. Taking the (fine-tuned) 
limit, we show that our Wigner functions do reduce to the expected ones 
on the half-line. This confirms that the Wigner transform should indeed 
be unmodified for this contact interaction. 



PACS: 03.65.-w, 03.65.Db, 03.65. Sq, 03.65.Nk 



1 



1 Introduction 



Point interactions and reflecting walls are known as contact interactions pQ. In 
quantum mechanics, they have been a subject of some interest lately-see [2] for 
some intriguing properties. Perhaps the simplest example is an infinite reflecting 
wall [J. 

Contact interactions are described by potentials with sharp features. A 
smooth interaction can be encoded in a potential feature of a certain width 
w. A contact interaction is obtained for zero width w, or, equivalently, infinite 
sharpness a := l/w. Physically, these sharp features should be understood in 
terms of the limit w — > (a — > oo). 

But can systems with sharp features be quantized after the limit is taken, 
or is it necessary to quantize before? Do we take a — > oo before or after 
quantization^ 

If one takes the sharp limit a — ¥ oo before quantization, one can rely on 
mathematical conditions to proceed. In operator quantum mechanics, one only 
needs to impose appropriate boundary conditions on wave functions in coordi- 
nate space. The boundary conditions conserve probability and can be under- 
stood as necessary for self-adjointness of Hermitian operators, like the Hamil- 
tonian, or extensions thereofo For the infinite reflecting wall, Robin bound- 
ary conditions are the only possibilities. They include the standard Dirichlet 
boundary conditions and the Neumann conditions as two extremal points in a 
one-parameter continuum of possibilities. 

It has been emphasized that the non-standard versions of such interactions 
should not be ignored, since they may describe interesting physics [6j [T] . How- 
ever, physical considerations such as symmetry (such as time-reversal invariance, 
e.g.) can eliminate possibilities in some cases [7]. Can the physical possibilities 
be restricted in other ways? 

Our point of view is that physically, zero-width (or sharp) features must be 
understood fundamentally as a — > oo limits of nonzero-width (smooth) ones. 
That is, the sharp case is an idealization, whose treatment should only provide 
a shortcut to the results obtained in the physical limit. 

1 One way to anticipate that it does make a difference is to realize that the classical limit 
h — > and the sharp limit a — ¥ oo do not commute. The wave phenomenon of non-Newtonian 
scattering [3] [4] makes that plain. For particle energy exceeding a discontinuous potential, 
there is a non-zero probability of reflection off the sharp feature, even though the process does 
not occur classically. Most strikingly, the probability is independent of Planck's constant, and 
so does not vanish as h — > 0. 

2 See [5] for nice expositions of the theory of self-adjoint extensions. 
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In that spirit, the infinite potential wall was described by a limit of the 
Liouvillc potential in [8] . The standard wall, with Dirichlet boundary conditions, 
was recovered. In this work, we will extend the result of [5] to the general case 
of Robin boundary conditions, by generalizing the Liouville potential to a Morse 
potential. 

In agreement with the results of others [T] , mass-dependent fine tuning is 
found to be necessary for non-standard walls to emerge. We believe that this 
fine tuning explains why standard quantum walls, with their Dirichlet bound- 
ary conditions, are standard. Non-standard walls are unlikely to be realized 
physically, because the required fine tuning is improbable^ 

Our original motivation came from phase-space quantum mechanics (a.k.a. 
deformation quantization) Q A further complication arises from sharp potential 
features in this context 13 H We therefore also examine Wigner functions for 
non-standard and standard walls, using the Morse potential. 

Dias and Prata [T3] treated the special (standard) case of Dirichlet boundary 
conditions for the Schrodinger wave functions. To describe the complication 
they found, let p(x,p) denote the Wigner function. At finite a, it can be found 
two ways. First, one can start from the wave functions ^{x), and build the 
corresponding density operator. Then a Wigner transform will yield the Wigner 
function; denote the result p a [ip](x,p). Alternatively, one can use the dynamical 
equations of phase-space quantum mechanics. The ^-eigenvalue equations can 
be solved, to yield p a [*](x,p). As long as a < oo, we must have 

Pa[^](x,p) = Pa[*](x,p) . (1) 



For the case a = oo, Dias and Prata found 

Poo[lp](x,p) ^ PooMfop) • (2) 

They then assumed that the Wigner transform p oo [-0](a;,p) was unaltered and 
added a boundary potential so that the ^-eigenvalue equations were compatible. 
That is, they modified the *-eigenvalue equations so that their solutions were 

Poo[*](x,p) = Poallflfap). 

In an effort to justify their somewhat ad hoc approach, alternatives to the 
^-eigenvalue equations were found in |14j . Dias and Prata then demonstrated 

3 This is reminiscent of the result of I10| , where the renormalization of a different singular 
interaction was shown to select a preferred self-adjoint extension. 

4 See for a review and 1 1 2| for a pedagogical introduction, e.g. 

5 Perhaps this is not surprising, since even the corresponding classical trajectories are 
continuous in configuration space but discontinuous in phase-space. 
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[T5] that the use of the alternative, so-called *-eigen-* value equations, had 
a certain equivalence to their treatment. Since the *-eigen-* value equations 
were derived, rather than postulated, those arguments provided an indirect 
justification of their procedure. 

More directly, the limit of the Liouville potential was studied in [5] . It was 
shown there that the Wigner transform of the wave functions with Dirichlet 
boundary conditions was indeed physical, as was assumed by Dias and Prata 
[13]. That is, 

lim p a [ip](x,p) = lim p a [*](x,p) = p^ [ip] (x, p) . (3) 

a— >oo a— >-oo 

The first equality was guaranteed, by (JlJ, but the second was not. If the limit 
had produced p^, [*] (a;, p) instead, for example, then the Wigner transform would 
have had to be modified, rather than the ^-eigenvalue equations. 

In this work, we will extend the result Q of |8] to the general case of Robin 
boundary conditions. 

Let us also mention that in [5], the connection was first made between self- 
adjoint extensions and the problem ([2]) found by Dias and Prata [13]. Subse- 
quently, those authors were able to show that the Hamiltonians that included 
the boundary potentials they introduced were indeed self-adjoint [16] . 

This paper is organized as follows. In the next section, we review the re- 
alization of Robin boundary conditions in limits of piece- wise flat potentials, 
following [9] Q]. There, mass-dependent fine tuning of the potential was found 
to be necessary to realize a non-standard wall, i.e. to avoid the standard Dirich- 
let boundary conditions. We point out that this fine tuning is equivalent to 
selecting a reflection resonance, as defined in [T7] , 

In section 3, the analogous calculation is carried out for a smooth Morse 
potential. The Robin boundary conditions are recovered, with the same kind of 
mass-dependent fine tuning already found in [9] [T] . We also show that reflection 
resonances are again selected in the smooth case. 

In section 4, Wigner functions for the Morse potential are considered. Using 
our solutions of the ^-eigenvalue equations, described in [TS], we show that the 
Wigner functions reduce to the expected ones [1] in the appropriate limit. That 
is, eqn. © is indeed obeyed. 

The final section is our conclusion. 
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2 Robin boundary conditions from a 
discontinuous potential 



Consider a non-relativistic quantum particle that is confined to the positive 
half-line with coordinate x, but is otherwise free. Its wave function must satisfy 
the so-called Robin boundary conditions 

^(0) + LV'(O) = (4) 

for some real length parameter L 6 (— oo,oo)U{oo}. The Robin, or mixed 
boundary conditions generalise the Dirichlet (L — 0) and Neumann (L — > ±oo) 
ones. They conserve probability and realize the self-adjoint extension of the 
Hermitian Hamiltonian H = p 2 /2m on the half-line. 

Though there is no mathematical reason other than simplicity to prefer 
them, Dirichlet boundary conditions are the most commonly applied. For that 
reason, infinite walls with other boundary conditions imposed are known as 
non-standard walls [T] . In this paper we investigate the physical motivation for 
standard and non-standard walls. 

The real wave function 

i>k{x) = sin(kx + (j)) (5) 

obeys the boundary condition if the phase is chosen so that 

kL — — tan</> . (6) 

It is appropriate for an unbound particle of energy h 2 k 2 /2m. For the same 
dynamics, one bound state also exists, with (unnormalized) wave function e~ x l L 
and energy —h 2 /2mL 2 , provided L > 0. 

The bound state provides the length scale L: its energy defines it, and its 
wave function has range L. This does not work for L < 0, however. A more 
democratic interpretation is provided by the Wigner time delay (advance) 

5t = 2% dE = -hk(l + k 2 L 2 ) ■ (?) 

for L > (L < 0)@ 

Let us now consider a particle moving on the (whole) real line with coordinate 
x and Hamiltonian 

H = p 2 /2m + V(x) . (8) 



6 See [T] and references therein. 
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A particle with energy 2mE = Ti 2 k 2 has a time-independent wave function ip(x) 
satisfying the stationary Schrodinger equation 

h 2 d 2 %jj(x) , . , . h 2 k 2 



2m dx 2 + V M*M =1^^)" 
We will show that Robin boundary conditions can arise from the limit of a 
smooth potential. This generalizes the derivation of Dirichlet boundary condi- 
tions from the sharp a — > oo limit of the Liouville potential V a (x) — 4^- e~ 2ax . 
In the context of deformation quantization, the latter result was obtained in [5] . 

To prepare for that calculation, we'll first study a discontinuous, piece-wise 
flat potential: 

!oo , x < , 

-££al{al + l), 0<x<l/a, (10) 
, x > 1/a . 

Here £, 1/a and 1/k are lengths, with k 2 > controlling the overall strength 
of the potential. Seba [9] showed that when this potential becomes an infinite 
wall as a — > oo, Robin boundary conditions are recovered. 

To see this, solve the Schrodinger equation piece-wise to get 

{0 , x < , 

sin (x^k 2 + k 2 at (a£ + 1)) , < x < 1/a , (11) 
A s'm(kx + (f>) , x > 1/a 

for an energy E = U 2 k 2 /2m > 0. Notice that the boundary conditions at x = 
are Dirichlet. Those at x = 1/a, however, are of the mixed type, i.e., Robin. 
We can therefore derive Robin boundary conditions at x = 0+ := lim Q _ ! . 0O 1/a. 
From the point of view of the physical wave function outside the resulting point 
interaction, it is the Robin (instead of the Dirichlet) boundary conditions that 
must be imposed. 

Matching the wave-function values and derivatives at x — 1/a, and taking 
the large a limit gives 



TT f 1 

l( n+ 2 



n G Z . (12) 



Then sin(K„^) = ( — 1)™, and we find 




if n _|_ 1)4 2k£ 
A = A n := + , tan0 = tan0„ = • (13) 

Comparing to ([5]), we get 



21 

ir 2 {~\) 2 



L = L n ■■= _ 9 ,_ I TT 2 (14) 
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for the Robin length scale. 

So the Robin boundary conditions arc found for x — 0+, but only barely: 
there are solutions only for a discrete set of values of re, indexed by the integer 
n. The strength of the potential needs to be finely tuned, tuned differently 
for different particle masses^ and the non-standard Robin boundary conditions 
arise for a very limited subset of possible parameters. 

What is the physical significance of the fine tuning? It selects a resonance. 
For this potential the probability of reflection is always one, but a reflection 
resonance can still be defined by a rapid change of tt in the phase shift [TT] . 
From the matching conditions we can derive 

tan(j/a) tan(fc/a + <f) 



(15) 

where j := [k 2 + n 2 a£(a£ + 1)] 1 ^ 2 . Demanding that = and selecting the 
maxima of leads to tan(j/a) = oo, or 



j/a k/a 

f-i-,of n 

Hep 



J 



[k 2 + K 2 al{al+l)] 1/2 _ ( 1 



- 7T 



neZ . (16) 



a a \ 2, 

In the a — > oo limit, the fine-tuning condition (|12p is recovered. 

Let us note that the reflection resonance condition (IT6|) corresponds to Neu- 
mann boundary conditions at x = 1/u, even before the a — ¥ oo limit is taken. 
Of course, the requirement (JT2]) for Robin boundary conditions does not select 
Neumann boundary conditions. Substituting (|12l) yields 

j - = ( n +^) T + (^n)" 1 + 0(a- 2 ) , (17) 

using (|14j) . This shows that the fine tuning is to near a reflection resonance; 
how it is approached in the a — > oo limit determines the Robin length scale L n 
and so the boundary condition that is realized. 

Let us also consider the bound states of the Seba potential in the a — > oo 
limit. For the negative energy states the wave function will decay exponentially 
in the interval x G (l/a,oo). Dividing the matching conditions for the wave 
function and its derivative yields 



2m\E\/hr + K 2 a£(a£ + 1) cot - d-2m\E\/tf + k 2 a£ (a£ + 1) 



= -\l2m\E\/Ti 2 . (18) 



7 For the standard wall with Dirichlet boundary conditions, this mass dependence is not 
present. 
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As before we can compare the coefficients in front of the different powers of 
a. The only possible energy is then E = — fi 2 K 4 l 2 /8m. Taking into account 
that (|12M14I) are needed for the Robin boundary conditions to arise, we obtain 
the correct energy —Ti 2 /2mL 2 in the a — » oo limit. The bound state energy is 
recovered from the Seba potential, for L > 0. We can also verify that, for the 
same values of the parameters, the (unnormalized) wave function of the unique 
bound state is e~ x / L in the limit. 

One criticism of these results could be that Dirichlet boundary conditions 
were assumed, not derived for the infinite wall (with no extra structure) at x — 
in the Seba potential. In addition, infinite potential walls are only idealisations 
of very high, but finite walls, and so the infinite wall should be treated as the 
limit of a finite wall. However, similar results were obtained later in [I] but with 
a finite wall, and no particular boundary conditions assumed. Robin boundary 
conditions were again obtained, with non-standard walls arising only when a 
mass-dependent fine tuning was imposed. We will therefore study here a non- 
sharp, or smoothed version of the Seba potential, rather than of the potential 
in [U. 

The authors of [T] speculate that a better choice than their piece- wise flat, 
discontinuous potential might eliminate the peculiar mass-dependent fine tuning 
required for non-standard walls. Presumably, it could also be argued to be pos- 
sible for the Seba potential [5]. We will find, however, that the mass-dependent 
fine tuning remains necessary in a smoothed version of Seba's potential. In 
retrospect, this should perhaps not be surprising, at least for Schrodinger quan- 
tum mechanics. The limit that squeezes and stretches the potentials into an 
infinite wall is so extreme, it seems unimportant whether the original potential 
has corners or is smoothed. 

3 Wave functions with Robin boundary 
conditions from a Morse potential 

To study how Robin boundary conditions arise as limits in deformation quan- 
tization, sharp potential features should be avoided. We will now carry out an 
analysis similar to that of the previous section, but for a smooth potential. The 
spectrum is first found for a potential with undetermined parameters. Then 
we consider a certain limit of the parameters, demanding that we recover the 
infinite wall, and that the states realized coincide with the eigenstates for the 
infinite wall, with Robin boundary conditions obeyed. 
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With its short range repulsion and longer range attraction, the smooth Morse 
potential 

V{x) = (e~ 2ax -be~ ax ) (19) 

is a rough approximation to Seba's. Besides a, two more parameters are needed 
- k determines the overall potential strength, b > the position of the well, 
and together they fix its depth. We will need to impose conditions on these 
coefficients in order to obtain Robin boundary conditions in the large a limit. 

The previous section indicates that we need only show that Robin boundary 
conditions apply at x = e, where epsilon is very small, but beyond the features of 
the Morse potential when a —¥ oo. For the unbound wave functions, therefore, 
we need only require that the relevant wave functions have the asymptotic form 
ip(x) ~ Asm(kx + <j>) as x —¥ oo, with <p variable. 

To do that we first need to find the unbound wave functions. Following 
Matsumoto [115] , we can solve the stationary Schrodinger equation for the Morse 
potential (|19|). The substitution tp{x) — <fi{z), z = exp(— ax), changes the 
Schrodinger equation into 



zV + Z(f>' + \ 



2mE 2 2, 2, 
s K Z + K Z 

h 2 



= . (20) 



This can be further transformed into canonical form (without a first derivative 
term) using the substitution <f>(z) = z" 1 ' 2 F{z). Changing the variables to 
y := 2nz/a leads to the so-called Whittaker equation, treated in [20], Chapter 
XVI: 



1 bn 1 1 

1 + o + S 

4 za y 



1 fik 
4 ~ \ ~a 



f = , (21) 



where f(y) :— F {ay /2k) and, as before, k — \/2rnE/h. The two linearly inde- 
pendent solutions are defined in [3T], p. 755. They are called Whittaker functions 
and can be expressed in terms of the Tricomi confluent hypergeometric function 
U{n, z) and the Kummcr confluent hypergeometric function M(/i, v, z):[§ 

M lm {z) = z ™+i/2 e -^/2 M(l/2 +m-l,l + 2m; z) , (22) 
Wi m {z) = z ™+i/2 e -z/2 U{l/2 + m-l,l + 2m; z) . (23) 



8 The Whittaker function Mi m {z) should not be confused with the Kummer function 
M(ji, v, z) in the above equation. Subscripts are used to denote the parameters of the Whit- 
taker functions in the literature, and the explicit bracket notation is used for confluent hyper- 
geometric functions. For further information involving the hypergeometric functions see |21j . 
p.753 and [22], p.503-506. 
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For our purposes, we only need the definitions of those functions 

Mfaw) = (24) 

Ufa, v- z) = r ^ ~ 1) z 1 - y M(l + m - ^ 2 - i/; 2) 

1 T ^~ v > Mfa,u;z) . (25) 



r(/*-i/ + i) 

Here we use the Pochhammer symbol fa) n := /x(/x + l)...(/i + n — 1), (/i)o := 1. 
Now the wave function can be written as 

1>k(x) = e ax ' 2 [c 1 M^^(y(x))+C 2 W^ My(x))] . (26) 

Imposing reality yields C\ — 0. The second term has physical asymptotic be- 
haviour: for large positive x it is sinusoidal with a phase depending on the 
potential parameters; for negative x far from the origin, there is the expected 
rapid exponential decay of a classically forbidden region. The wave function is 
therefore 

i>k(x) = Ce ax/2 W^_^ (^ e ~ aX ) • ( 27 ) 

With the help of equation (|25l) we can rewrite this result in a form similar 
to that given by Matsumoto in [TO] for a Morse potential with b — 2. The wave 
function is manifestly real in this form: 

4>(y) = Ce-y^Ay^M(\- h ^ + ^,l+ 2 ^;y 



2a a 

+ Ce-y/ 2 A* y-*/°M (i - - l l,l-™i.y) , (28) 



with C a real normalization constant, and 

p/ 2ik \ 

A = v( i_^_ m ■ (29) 

V 2 2a a) 

Let us now examine the asymptotic behaviour of the wave function and how 
it depends on the parameters. In the limit a — > 00, exp(— y/2) ~ exp(— e~ ax ) 
approaches the step function, so the dynamics will be restricted to the positive 
half- line. The limit x — > 00 corresponds to y — > 0. Using (|28p and (j2"4")l we 
obtain 

ip(x) ~ C|i| cos \kx-axg(A)] . (30) 



10 



The phase can be calculated from and Euler's infinite product formula 
1 



»7« 



T(u) 

A short calculation shows that 

7T 7& 



n[ 



i 



&rg{A) 



a 



(31) 



(32) 



OO 

E 

n=0 



a(n + 1) 



tan 1 


2k 


+ tan" 1 


a(n + 1)_ 







2A- 



(2n+ l)o - 6* 



Apart from the tt/2, all terms will vanish in the a — > oo limit, except those of 
the form tan -1 [2k/ ((2n + l)a — 6k)]. For one such term to survive the limit, 
we need k = 0(a 1 ). If the strength k does not have this form, we will recover 
Dirichlet boundary conditions, i.e. the standard wall. Now, since it is bn that is 
relevant, we let b absorb the proportionality constant, and use k = a + 0{a ). 
Finally, because the terms of order a and lower will not affect the results, we 
drop them, and put n = a from now on. 

In order to realize Robin boundary conditions (|4|), the parameter b must be 
of the special form 



b= (2n + l)- 2L~ l /a + 0(a~ 2 ) 
Here L is a fixed length, independent of a. Then we find 

kL = tanarg(A) 



(33) 



(34) 



in the large a limit, so that the wave function ([28]) does indeed satisfy the Robin 
boundary conditions Q. 

At large a the term 2_L _1 /a is negligible compared to the other two. While 
the parameter b approaches an odd integer the second infinitesimal term is cru- 
cial. Apparently, we need to fine-tune the parameter b to recover the Robin 
boundary conditions. A version of this phenomenon has already been encoun- 
tered in |Tj where the parameters can only take very limited values. The authors 
argue that fine tuning may be a result of the particular choice of potential they 
are using, possibly because it is not smooth. Since our analysis, using a smooth 
potential, produces a version of fine tuning as well, fine tuning cannot be related 
to discontinuity alone. 

For non-standard walls, we must fine-tune the parameters so that we are 
near Neumann boundary conditions. Notice that this is precisely as it was for 
the Seba potential of sect. 2 (see (|17l) and nearby). Put another way, the 
fine-tuning is again to a reflection resonance, or slightly off its peak. 
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Let us now consider the bound states. Their wave functions are given in [23] 

as 

i>(x) oc exp(- Ke - Q " /a)e- a{ "- bK/2a+1/2)x L b l , K/a ^- 1 (2 K e- ax /a) , (35) 
where 

*<•>-£<-'>-(;:*)£ (36) 

m=0 v 7 

are the associated Laguerre polynomials L^(x). The energies are 

h 2 2 

E v = --^(v-bK/2a+l/2) 2 , (37) 
2m 

for integer v s [0, [6/c/2aJ], where [aj is the smallest integer less than a. 

Consider now the a — > oo limit. The Laguerre polynomials are normal- 
ized to one at zero argument, and exp(— ne~ ax /a) turns into the step func- 
tion. The only term that remains to be analyzed is e - Q ( iy - bK / 2Q + 1 / 2 ) a: . Clearly, 
—a{y — bn/2a + 1/2) must be a negative constant (independent of a) so that 
we have a normalizable wave function that does not disappear in the large a 
limit. Again we can set k = a, and the solving for b yields precisely equation 
(|33[) . Analyzing the bound states provides an alternative way of deriving the 
fine tuning condition. 

On the other hand, let us give particular values to the constants in (I33[) . 
i.e. fix b. All the bound states will vanish for x > except the one that has 
highest quantum number v = \b/2\. This is because the maximal integer will 
cancel the integer part of b and leave only the fine tuning part —2 (La) -1 . The 
bound state wave function ~ e~ x l L will be recovered with the correct energy 
-h 2 /2mL 2 . 

To summarize, in this section we demonstrated that the a — > oo limit of the 
Morse potential (|19p can be used to generate Robin boundary conditions. Fine 
tuning is necessary, however: the parameter b must be an odd non-negative 
integer plus a term with 1/a asymptotics that determine the length scale L 
of the Robin boundary condition. If the fine tuning is absent or if the integer 
part of b is not an odd integer, then we can only recover Dirichlet boundary 
conditions, i.e. standard walls. A new observation is that the fine-tuning selects 
a reflection resonance. 

Notice that the definition of k involves the particle mass. So, if the particle 
mass changes, so must the potential. The fine tuning required is also mass 
dependent. 
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The situation is very similar to that for the discontinuous Seba potential [5] 
treated in sect. 2, and to the results of HQ. The mass-dependent fine-tuning 
that is necessary for non-standard (Robin boundary condition) walls seems to 
be more than an artifact of the choice of potential. In particular, just smoothing 
out the discontinuities of a piece-wise flat potential is not sufficient to avoid this 
property. As already stated, this is perhaps reasonable in hindsight: it seems 
that the limit that squeezes and stretches the potentials into an infinite wall is 
so extreme that it is unimportant whether the original potential has corners or 
is smoothed. 



4 Wigner functions and Robin boundary 
conditions with a Morse potential 

In phase-space quantum mechanics (see JTJ [T2], e.g.), the Wigner function 
p(x,p;t) encodes all measurable information about the quantum state of a sys- 
tem. It satisfies the equation of motion 

ihd t p(x,p,t) = [H,p(x,p,t)]< , (38) 

where [H, p]^ = H * p — p* H, and the Moyal ^-product is defined by 

* = cxp | y ( d x d p - d p d x ^ | . (39) 

It can be expressed as a linear combination of stationary Wigner functions with 
time-dependent coefficients: 

p(x,p,t) = £ C ELER e-^- E «Wp ELER (x,p) . (40) 
Here „ denotes the Hamiltonian *-eigenfunction that can be found by solv- 

' E L E R b J 

ing the system of equations: 

H *P El e r ^P^ = El Pe l eS x 'P) ' ( 41 ) 
P ElEr (*,p) *H = E r P ElEr (x, P ) . (42) 

Alternatively, the Wigner transform 

/oo 
dy e w (x + hy/2\E L ) (E R \x - hy/2) (43) 
-oo 

allows them to be determined from the wave functions, if known. For smooth 
potentials, the resulting Wigner functions are known to agree. 
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For discontinuous potentials, however, that is not necessarily the case [T5] . 
For the infinite wall (or a particle confined to the half-line), Dias and Prata 
[T5] showed that the Wigner transform of the density operator only satisfies 
the ^-eigenvalue equations if the free Hamiltonian is modified. No independent 
motivation was given for the change to the Hamiltonian, however. It was also 
assumed that the Wigner transform itself did not need to be adjusted. 

An independent motivation was first suggested in [5]: self-adjointness of 
the Hamiltonian. The free Hamiltonian on the half-line is not self-adjoint. It 
does have self-adjoint extensions, however, and these correspond precisely to 
the possible boundary conditions (g]) (see [5], e.g.). Subsequently, the self- 
adjointness of the Dias-Prata modified Hamiltonian was demonstrated in |16j . 

Here we are concerned with the assumption of an unmodified Wigner trans- 
form. That is, does the unmodified Wigner transform of the density operator 
provide the physical Wigner function? In jS], we answered in the affirmative, by 
treating the infinite wall as the limit of a smooth, Liouville potential. Only the 
standard Dirichlet boundary conditions were recovered, however. Here we will 
show that non-standard walls can be realized in a similar way, using the Morse 
potential, and that the naive Wigner transform does indeed work, for all Robin 
boundary conditions, describing both non-standard and standard walls. 

The Wigner transforms of the density operator elements relevant to Robin 
boundary conditions have already been computed, in [¥]. For x > 0, using the 
wave functions (pj [5]), we find: 

sin [2(p/K - k)x] sin [2(p/h + k)x] 



Poc[^]{x,p) OC 



(p/h - k) ( P /h + k) 



t- x sin(2a; p/h) . , 
2cos(2fcz;-4) — K Z . 44 ) 

p/h 



In addition, the ^-eigenvalue equations (|4T1 1421) for the Morse Hamiltonian 
have also been solved in [18] . We must take their limit a — > oo as described in 
the last section, and compare with (l4"4"l) . 

In [18] , the ^-eigenvalue equations (|4il |42)) were solved for the Morse poten- 
tial using a Mellin transform and factorization^] Writing 

h 2 k 2 h 2 h 2 

the result was of the form 

9 The method used there should be useful for any potential that is a polynomial in an 
exponential. 
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f c+zoo 



The factors can be written using 

At+iki/2afn(_ 2ikr/a) 

Mt, **) oc r(1/2 . 6/ L, fcJ y) r( - 2t + * w x (47) 

2^1 (1/2 - 6/2 + ifa/a, -2t + iki/a; 1 + 2 ifa/a; 2) + 
4*-ifci/2ar(2ifc//a) 



-r(-2t - ifcj/a) x 



r(l/2 - 6/2 + tfc//a) 

2-F1 (1/2 - 6/2 - ifcj/a, -2* - ifc//a; 1 - 2ifc//a; 2) , 

with I — L. R. Here we have defined ^4 = 4 J ' c - r / Q Ji, with A as in (|29p. and used 
the identity T(2z) oc 2 2z - 1 T(z)T(z + 1/2). This is the general solution of the 
^-eigenvalue equations of phase-space quantum mechanics, for unbound states 
in a Morse potential (fT9]) with arbitrary real parameter 6. 

Following [5], we use the residue theorem to find the limit of the Wigner 
function when a — > oc. Since the calculation is straightforward but lengthy, we 
omit the details. The integrand of (|46p has 4 terms, one proportional to A 2 , 

one to A* 2 and two to \A\ 2 . The |A| 2 -terms yield contributions proportional to 

[ e 2«(p/?i-fc) _ e -2ix(p/n-k)y^ n _ and [ e 2ix( P /h+k) _ ^ixtp/h+Qyfjp/fi + k y 

The A 2 -term and the A* 2 -term yield 

^^2i A+2ixp /U — 2ixk ^2i arg A — 2ixp/h — 2ixk~\ I ^ 



and 



g — 2i arg A J r2ixp / U+2ixk ^ — 2i arg A— 2ixp/h J r2ixk j 



Note that all the terms arising from residues at i(±p/7i±' k)/2a + 1/2 produce 
decaying exponential factors and therefore do not contribute. Also, the Gauss 
hypergeometric function 2-F1 (a, 6; c; z) is analytic with respect to its second ar- 
gument and 00 is its only singularity. The contributions from 2-Fi(a, 6; c; z) will 
manifest themselves as a multiplication by constants in all cases. In particular, 
for those terms that survive in the limit of interest, the constant is 1. 

The algebra can now be completed to reproduce the Wigner function for 
an infinite, but possibly non-standard wall as we hoped. The Robin boundary 
conditions are indeed recovered using the Morse potential. We have outlined 
how the calculation is done for diagonal elements of the symbol of the density 
operator, but the non-diagonal case works in similar fashion. 
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5 Conclusion 



Let us summarize our results. 

In section 2, we reviewed Seba's analysis [3] showing that Robin boundary 
conditions (for wave functions) could be realized by a limit (a — > oo) of a dis- 
continuous, piece- wise fiat potential, eqn. (110[) . We pointed out that standard 
walls (Dirichlet boundary conditions) are generically realized in the sharp limit, 
and non-standard walls arise only if a mass-dependent fine-tuning (|12p is im- 
posed; these observations are in agreement with those made in pQ , for the limit 
of a similar, but everywhere finite, potential^ We observe that the parame- 
ters are fine-tuned to a reflection resonance in the limit. If the fine-tuning is 
imposed, then non-standard walls can be realized, and the Robin length scale L 
is determined by exactly how the limit resonance is approached (see eqn. (jTTJ)). 

The analysis of the piece-wise flat Seba potential was repeated with a qual- 
itatively similar, but smooth potential, the Morse potential of eqn. (fT^j) . Re- 
markably, the results were almost unchanged. Analysis of both the unbound 
and bound states yielded a mass-dependent fine tuning ([55]) required for non- 
standard boundary conditions. Again, a reflection resonance is selected by the 
fine tuning, and how the resonance is approached in the a — > oo limit determines 
the precise boundary conditions realized, i.e., the Robin length scale L. 

As mentioned in the introduction, the infinite reflecting wall is perhaps the 
simplest example of a so-called contact interaction. For such, the interaction im- 
poses boundary or matching conditions, such as the Robin boundary conditions 
on the half line. Alternatively, the same conditions can be found by demanding 
that the Hamiltonian or its extension be self-adjoint (see [5])- 

Here we have assumed that contact interactions can only be realized phys- 
ically as limits of smoother, less localized interactions In agreement with 
[5J [T] , mass-dependent fine tuning was found to be necessary for non-standard 
walls to emerge. We therefore believe this explains why standard quantum 
walls, with their Dirichlet boundary conditions, are standard. Non-standard 
walls are unlikely to be realized physically, because the required fine tuning is 
highly improbable. 

It would be interesting to see if the realizations of other contact interactions 
as limits require similar fine tuning, and if there are other so-called standard 
boundary/matching conditions selected that way. 

10 The independent derivation of the fine-tuning condition from a study of the bound states 
(sec eqn. (JTSJ), instead of just the continuum, is perhaps new. 

11 Clearly, we do not consider effective motion for a radial coordinate r G [0, oo). 
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Finally, our primary motivation came from phase-space quantum mechanics, 
or deformation quantization. In that context, we made some progress on solving 
the dynamical equations of Wigner functions, reported in [18 . Here, we were 
able to demonstrate that in the sharp limit, our Wigner functions become those 
constructed by the Wigner transform from wave functions with Robin bound- 
ary conditions. Specifically, we showed that when a — > oo, equations (prol ET71) 
reduce to the expected Wigner function, eqn. [J, x > 0. This justifies 

the assumption that the Wigner transform is unmodified for these examples of 
contact interactions. 
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